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We study some generic aspects of the winding angle distribution around a point in 
two dimensions for Brownian and self avoiding walks (SAW) using corner transfer matrix 



p • and conformal field theory. 



1. Introduction 

The topological constraints on Brownian walks are a long standing subject of interest 
in mathematics and polymer physics. A particularly fruitful example is the study of the 
winding angle 9, the continuous angle swept by the movement around a set of points 
(curves) in two (three) dimensions i. In a plane, the probability distribution at large time 
or length I for the winding angle of a Brownian walk around any point O was determined 
by Spitzer [^j as 
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^ Since all results are even in 6 in what follows we generally consider 6 > 0. Otherwise, 
substitute \6\ for 9 



In what follows the length of the walk is usually thought of as the number of steps in 
a discretized version. Because of the logarithmic dependence on / the various possible 
definitions of all proportional to each other, are equivalent. This of course is compatible 
with angles being dimensionless. 

The law ( p. . 1|) is substantially modified when a cutoff is introduced so that the walk 
cannot get arbitrarily close to O. One finds then 
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The standards methods of approach to these problems are based on diffusion equations, 
functional integrals or refined probability theory . 

More recently the same problem has been considered for self avoiding walks (SAW). By 
applying methods derived from Coulomb Gas representations and conformal field theory 



, the equivalent of (|1.1| ) or ( p..2|) (for SAW there is no difference between the two cases as 
the walk provides a natural UV cutoff, see figure 1) was determined |^] 
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The purpose of this note is to study the relation between (|1.2|) , ( p..l|) and ( [17^ ) 
further, and to compute a related distribution for SAW that will also depend on the 
variable 6 / log/. The analysis uses corner transfer matrix ideas and conformal field theory. 
It can be considered follow up of H . 

Indeed, the most striking difference between the Brownian and self avoiding cases 
lies in the variable that is either 6 / \ogl or ^/(log/)^/^. The intuitive explanation of this 
difference is that the distribution in the self avoiding case is mainly determined by the 
excluded volume of the walk already wound around O (figure 1), while in the brownian 
case, there is no such effect, and the distribution is determined fully by the entropy loss that 
arises from the constraint of winding . To suppress this major difference and concentrate 
on curvature related entropy we simply put the SAW on the multi sheeted Riemann surface 
for the function In z, discretized if necessary. For a walk of length / we define the probability 
of winding angle 6 by the relative number of configurations that sweep a total angle 6, 
and for |^| > 27r have end points on different sheets. It is important to realize that this 
problem would be identical to the usual winding angle problem for the brownian case 
(since there is no interaction between different parts of the walk one can just collapse the 
staircase onto the plane). 

To determine the probability distribution we can proceed in two ways. 



2. Conformal mapping 

1. First we consider a lattice SAW on a strip of length L and width W. The boundary 
conditions are periodic (free) in the L (W) direction (figure 2). On an infinite lattice the 
number of SAW of length 1 grows as oc /i', where fx is the lattice effective coordination 
number. On the strip, the number of configurations for SAW winding around the periodic 
strip is denoted by We then consider the generating function 

Zl=5]o^%-^ (2.1) 
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In the limit where L/a,W/a — > oo, L/W remaining finite (a is the lattice spacing), Z is 
determined by conformal invariance and Coulomb gas mappings arguments. It is important 
to notice that there is no bulk term, the reason being that SAW have a fractal dimension 
less that two. The rest of the procedure is standard, but we recall it very briefiy for 
completness. First, we choose for technical reasons the honeycomb lattice, and reformulate 
the polymer problem as the limit n ^ of an 0{n) model. Partition and correlation 
functions of this 0{n) model are expressed as sums over self avoiding mutually avoiding 
sets of loops and open walks, with a weight (3 (the inverse temperature) per monomer and 
n per loop 0. For n = 0, the critical temperature is (3c = iJi~^ ■ Next, one reformulates the 
0{n) model in terms of the Izergin Korepin (IK) vertex model 0, the correspondence 
between the two being conveniently carried out using quantum group symmetries 



The continuum limit of the polymer problem is then worked out using a combination of 
sectors of the continuum limit of the IK model, which is described by a Gaussian model 
with charge at infinity [0. The net result for the lattice model is the following. The 
partition function where all non contractible loops on the cylinder have weight n = has 
the expression 
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Zo = - dplCp + 'ip+i/2A^p+i/2^ (2.2) 

where /Cg is the partition function of the IK model in the spin s (integer or half integer) 
sector, and dg is the quantum dimension 
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with 

n = -t-t~^. (2.4) 



At t = z, using the additional level coincidences induced by Utsl{2) representation theory 
at g a root of unity |jT^ , one checks that Zq = 1 (see fl^). The computation of Zi follows 
from the correspondence between the IK model and polymers as 
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In the continuum limit (L/a, W/a oo, L/W finite) we use the standard result that [jT3 
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and P{q) = Y[n>o(^ ~ 1^)^ ^r,s denotes the standard conformal weights in the Kac 
notations [jl^, with, in the polymer case, h^^i = ^^^~24 ~^ ■ We obtain then 
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One can also consider the same quantity Z2 defined for two SAW, mutually avoiding, that 
wrap around the cylinder. One has 
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2. In the limit L, ^ 00 the lattice strip can be considered as continuous, and Z 
is a partition function for the conformal field theory that describes polymers. The various 
results regarding this conformal field theory can be followed through conformal mappings. 
We consider the mapping z' ja = exp(z/a) where z is the complex variable in the strip 
plane, the time direction L being along the purely imaginary axis. We then obtain the 
same generating function but for SAW wrapping on a staircase (with the two end "lips" 
identified, see figure 3) of inner radius r and outer radius R in ratio W ja = log(i?/r), and 
angle 9 = L/a: 
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To show this more precisely it is necessary to identify Z with the trace of a transfer matrix, 
and then use known results for eigenvalues, whose properties in the mapping are deduced 



from the ones of correlation functions of the theory |15] . We explain this in more details 
in the next section. 

It is important to realize the exact meaning of the angle 9: our results make sense when 
r/a, R/a oo, R/r finite and 9 —>■ oo, 9/W finite, since on the strip we had L/a —>■ oo. 
The above result still holds when r remains of the order of the lattice spacing, while R, 9 
still are very large []T^ |]T^[]T7|. This is due to the logarithmic dependence of W on the 



radiuses ratio: we discuss a similar case in the next paragraph. In the following we take 
r — a. 

3. Third we must discuss what becomes of this generating function when evaluated 
in the scaling region /3 ^ (in fact we will exclusively deal with the high temperature 
part /3 < The behaviour for low temperatures is quite different because of symmetry 
breaking fl^). Once again we can rely on the 0{n) model analysis. Recall that the limit 
/3 — > is the approach to the critical point, where the correlation length diverges as 

^oc(^-l-/3)-^ (2.11) 
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where v is also the exponent characterizing the size of the SAW (< {Rg/clY >oc P'^ ^ Rq 
the radius of gyration). 

Now, in most problems, the knowledge of some quantity at the critical point is not 
enough to determine this quantity away from it. Usually the functional dependence is 
constrained by scaling arguments, but a non trivial function has still to be determined. The 
present problem turns out to be simpler. Indeed the usual finite size scaling hypothesis, 
valid in the limit R/a,^/a oo, involves replacing for each of the terms in the above 
sums R by Rf{^/R) where the function / depends on the term considered. One needs 
/(oo) = 1 and also, to suppress the R dependence in case the physics is dominated by the 
much smaller scale ^, f{x) ~ /qx, x — |]T^. Now suppose we consider precisely the latter 
regime. Since ^ still has to be very large compared to the lattice spacing and the dependence 
is logarithmic, /o disappears from ( p.8| ) ( p.9|) as log(/o^/a) ~ log(^/a), ^/a ^ oo. Therefore 
we conclude that in the scaling region, (|2.8|) ( p^ ) still hold with R/a replaced by ^/a: 
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In effect, the IR cutoff at criticality is exactly replaced by the correlation length in the 
scaling region. 



4. So far we have discussed generating functions only, while results of most statistical 
interest rather deal with SAW of large but fixed length. The standard procedure to get 
such results involves inverse Laplace transform |]T^ , and is usually quite complicated for 
full probability distributions. Here again, the situation is simplified by the logarithmic 
dependence. Call the number of SAW of length / on the staircase. The scaling function 
F we are looking for is defined by 

^ ^,H-\lnlrF (^^y (2.13) 

The form of ( p.l3|) is explained as follows. For 9 small we should recover the usual result 
for a self avoiding loop, taking also into account the combinatorial fact that our loop can 



assume any position provided it encircles the origin . This determines that the power 
in O^. The exponent a is not known; it does not affect our results, and we set it to 
zero in the following. The generating function we are computing is then 

Zf- = J]0f/3' - / ^ (i^) exp [/(/3/. - 1)] J. (2.14) 

Since in the region of integration / >> 1, F in ( |2.14D varies very slowly over the domain 
where the other term in the integrand takes appreciable values. This domain is around 
^ ^ i-p^ ■ Since F depends logarithmically on I and we consider the behaviour as — fx~^, 
we conclude that up to proportionnality constants zf'''^"'^ oc F — j^^^zrg^ ie the fixed 
length distribution F is obtained by replacing in zf^""^^ ( (p.8| ) or (|2.9| ) ) ln(^/a) by In P. 
In other words, the IR cutoff is now given by the average linear size of the chain. We find 
thus 

^ (i^) " £ n(-l)"<,(-')'--)/«, (2.16) 
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where 

q = exp(-^). (2.16) 
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The distribution function F has the following asymptotic behaviours 

F{x) oc exp(— 57rx/8), x ^ oo; F oc — , x ^ 0. (2-17) 
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where the last result follows from Poisson resummation. The divergence at small argument 
arises from the fact that in our computation the origin of the SAW is not fixed, introducing 
an additional degree of freedom. See section 4. 



3. Corner transfer matrix 



We can also derive this distribution using corner transfer matrix ideas, which are a 
httle hidden above. 

We now work in the plane from the start. We can still use the polymers, 0{n), IK 
model correspondence. The main difference with the computation on the strip is that the 
evolution is described by a corner transfer matrix (CTM) instead of a row to row one. It 
is important to notice that the CTM has also the property of commuting with Utsl{2): 
this is most clearly seen in the hamiltonian limit where for the row to row transfer matrix. 



H = ^H„ and for the CTM, K = ^iH, ^ where [H„Utsl{2)] = 0. As a result, the 



generating functions Z^^""^^ have expressions similar to the ones on the cylinder, like eg 

zf''- = J2{p+i)(-rc,+^/2, (3.1) 

where now Cg is the partition function of the IK model in the "radial spin s sector" 



2^, £ = Tr (tctm) ■ Similar expressions appear in the computation of spontaneous 



magnetization for the Q > 4 state Potts model [|2| . In corner transfer matrix calculations 
9 is usually a multiple of tt/2, but in the large angle limit we consider this becomes 
irrelevant. 

Now the problem is to evaluate directly (without using the chain of arguments of the 
previous section) the continuum limit of C when (3 fJ'~^,d oo. When ^/a = oo it is 
known that the IK model corresponds to a free bosonic theory with action 

^ = /-£^^^^'' ^'-'^ 

where for polymers g = 3/2, (p is the usual free field variable, dual to the arrows of the 
IK model in its solid on solid reformulation with normalization Sep = 2n x spin. Now 
when ^/a >> 1 this action has to be modified by adding suitable terms that make the 
model massive. Suppose first we were dealing with a free massive boson, ie SA = w? J 0^, 
m = 1/^. Suppose we consider a sector of radial spin s ie with boundary conditions 

(/)(r = oo) — (^(r = a) = —2'its. (3.3) 

The classical solution (ie satisfying (A + m'^)(f)ci = 0) to this equation is 
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where Kq, Ki are the usual Bessel functions. One has 
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We now have to evaluate the contributions of these two terms to the action, ie compute 
9 J rdr in the limit where ma << 1. Using 



Ko{x) oc ln(x), Ki{x) oc x ^ 0, 
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it is easy to see that the dominant contribution comes from the kinetic term, behaving as 
l/ln(ma) while the mass term gives a contribution of order 1/ In^ (ma). The classical part 
therefore behaves as in the analysis of the critical model on the periodic strip, with the 
now familiar substitution L/W —>■ 9/ln{l/ma). Therefore we conclude that in the limit 
ma —>■ 0, 9 —>■ oo, 9/ In(ma) finite 



(3.^ 



Now we must determine Zq^- This is most easily done by studying the diffusion equation 
for Brownian motion on our Riemann surface (not surprisingly, an identical problem occurs 
in the study of the winding angle distribution for Brownian walks [H ). Calling G the Green 
function one finds (recall that the region of radius a around the origin is excluded in our 
problem) 
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dv (3.9) 



Using contour integration, this can be rewritten as an infinite discrete sum over the zeroes 
(in variable v) of the denominator, corresponding to a sum over eigenvalues of tctm- In 
the limit ma —>■ these zeroes occur at 
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contributing to the angular dependence of the Green function by a term exp[7r^n/ In(ma)]. 
We have therefore determined the eigenvalues of the CTM, up to degeneracies. We now 



can rely on a mode analysis and the study of the ground state energy in to conclude 
that 

-1/24 

= (3.11) 

where q = exp[— tt^/ ln(l /ma)]. Although the correct action describing polymers in the 
scaling region is rather the sine-Gordon model, it is easy to see that in the limit we are 
interested in, this does not make any difference: all that matters is the kinetic term and 
the presence of a length scale m = 1/^. Therefore combining ( p.8|) and ( |3.11J ) we recover 
the results of the previous section. 

The careful reader may notice that in usual corner transfer matrix computations there 
are additional terms, usually made of elliptic functions, that multiply the exponential of 
the corner transfer matrix hamiltonian . In fact they disappear in the scaling limit 
we study here. Finally, computations on an integrable version of the off critical polymer 



problem have recently appeared, which confirm our results |23 



4. Comments 



1. The distribution ( |1.2|) holds for a problem where the origin and end points of 
the Brownian walk need not be on top of each other. It looks quite difficult to obtain a 
similar quantity for SAW; instead of taking a trace of the transfer matrix, one also has 
to insert matrix elements of boundary operators, which are known in principle but quite 
complicated. It is easier to derive the winding angle distribution for Brownian loops. One 
finds 
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which diverges at small argument, like (|2.17|) did. 



(4.1) 



2. The distributions ( [4.1|) and (|1.2| ) for Brownian walks and (|2.15|) for SAW are similar 
in nature. To see this more explicitely expand ( |1 . 2| ) as 



X 



29 



7rexp(— Tra;) j ^^(— )"" exp(— yrna;) j = tcqb j ^^(" 



(4.2) 



^n=0 



Kn=0 



where 
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an expression similar to ( |2.12|) since the exponent v for Brownian walks is equal to 1/2. 

The UV cutoff necessary to get (|1.2|) and not ( |1 . 1|) appears now very natural from the 
corner transfer matrix or conformal mapping point of view. Distribution ( |1 . 1| ) arises when 
this cutoff disappears and the CTM spectrum is not quantized any more: this situation is 
not met in most other lattice models (like the 8 vertex model) where the spectrum of the 
CTM is discrete, and in the continuum there is always an implicit UV cutoff like in the 
SAW case. 

Let us emphasize that the integer spaced structure of the CTM spectrum away from 
the scaling region in integrable models has nothing to do with the above arguments. 



We also do not believe that this integer spaced structure is related to quantization on 



angular momentum as claimed in : this is especially clear in the scaling region. 



3. The same procedure could also be applied to revover the result of ||^] concerning 
a SAW on the plane that makes many turns, building in that fashion a core of forbidden 
region aroud O (figure 1). In this latter case one deals indeed with the plane, not the 
infinite staircase so LjW — and the winding angle appears in the spin variable s. Since 
the dimension of the ground state of the spin s sector grows like s^, this explains why the 
distribution ( |T73| ) depends on 6'^/lnL Let us take this opportunity to complete a point 
that was overlooked in . The scaling form of a row to row transfer matrix eigenvalue 
as deduced from conformal invariance A ^ exp(— tt/i/I^) holds only in the limit — > oo 
with s finite. To establish the Gaussian distribution of one needs to apply the same 
formula for s ^ oo too, but s oc ^JW so s << still. We do not know any proof that 
this is legitimate, but we can give arguments based on the correction to scaling analysis 



2^ . Considering for instance the analytic corrections induced by the terms L-2L-2 and 



L^27 -^-2 (where L_„ denotes the usual Virasoro algebra generators [|4| ) one finds that 



the scaled log of eigenvalue behaves as + so even as s oc \/W the result still holds. 



4. In conclusion, the winding angle distributions of the Brownian, SAW and probably 
other geometrical problems (like percolation) around a single point can be deduced from 
conformal field theory only, and are given by some kind of theta function. Notice that 
such distributions make reasonable sense as statistical problems only because the bulk free 
energies vanish. Otherwise the results would always be dominated by very large angles. We 
do not really believe that our distribution ( |2.15|) on a staircase has much practical interest. 
However it is amusing to notice that it can be reformulated slightly differently. "Unfold" 



the discretized periodic staircase on the plane as in figure 4. Locally this is like the plane, 
but there is curvature at the origin. Since in our SAW problem all edges are treated as 
having the same length, we consider this lattice as the triangulation of a manifold with 
curvature located at the vertices, pi = tt ^"'^' where qi is the number of neighbours of 
vertex i (with the surface element Si = ^ this satisfies Gauss Bonnet theorem for closed 
manifolds). Distribution ( |2.15| ) appears thus as the distribution for a SAW "trapped" 
around a point of strong negative curvature, in the limit where /, p oo. Maybe in that 
form will it prove more useful. 
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Figure Captions 



Figure 1: A SAW on the plane winding around the origin. 

Figure 2: Periodic strip. 

Figure 3: By conformal mapping the strip maps onto a periodic staircase. 

Figure 4: A discretized staircase "unfolded" on the plane with negative curvature located 

at the origin, and a self avoiding loop that winds around it (with 9 = 8^ = An). 
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